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Magnetism in the three-dimensional layered Lieb lattice: 
Enhanced transition temperature via flat-band and Van Hove singularities 

Kazuto NodaQ Kensuke Inaba, and Makoto Yamashita 
NTT Basic Research Laboratories, NTT Corporation, Atsugi 243-0198, Japan 

We describe the enhanced magnetic transition temperatnres Tc of two-component fermions in 
three-dimensional layered Lieb lattices, which are created in cold atom experiments. We determine 
the phase diagram at half-filling using the dynamical mean-field theory. The dominant mechanism of 
enhanced Tc gradnally changes from the (delta-functional) flat-band to the (logarithmic) Van Hove 
singularity as the interlayer hopping increases. We elucidate that the interaction induces an effective 
flat-band singularity from a dispersive fiat (or narrow) band. We offer a general analytical framework 
for investigating the singularity effects, where a singularity is treated as one parameter in the density 
of states. This framework provides a unified description of the singularity-induced phase transitions, 
such as magnetism and superconductivity, where the weight of the singularity characterizes physical 
quantities. This treatment of the flat-band provides the transition temperature and magnetization 
as a universal form (i.e., including the Lambert function). We also elucidate a specific feature of 
the magnetic crossover in magnetization at finite temperatures. 

PACS numbers: 67.85.-d, 71.10.Fd, 71.27.-l-a, 75.10.-b 


I. INTRODUCTION 

Phase transitions, such as magnetism and supercon¬ 
ductivity, are of fundamental interest in lattice fermions. 
As a common feature, transition temperatures Tc are usu¬ 
ally given as a function of the density of states (DOS) at 
the Fermi energy p{Ep) and an interaction U, TcjW oc 
^-t/p(Ep)u X, where bandwidth IF is a unit 

of energy. A singularity located on the Fermi energy 
[p{Ef) —>■ oo] changes this functional form, which could 
greatly increase the Tc- For instance, the logarithmic Van 
Hove singularity (VHS) induces a characteristic depen¬ 
dence: Tc/W oc [l|. Recent studies proposed 

the emergence of another interesting delta-functional sin¬ 
gularity, which we call the flat-band singularity (FBS), 
at the surface of a layered graphene 0 or of a topologi¬ 
cal material Q. This FBS is also expected as the origin 
of the higher transition temperature: Tc oc U. Although 
these singularities have attracted attention, we still lack a 
comprehensive understanding of these singularity effects 
on phase transitions. 

Cold atoms in an optical lattice, where we can control 
lattice geometry and resulting DOS singularity an , pro¬ 
vide opportunities for studying the singularity effects in 
bulk systems. In particular, successful creations of two- 
dimensional (2D) optical lattices with singular DOSs, the 
Kagome Q and Lieb (line-centered-square) [3] lattices, 
have activated theoretical studies on phenomena related 
to the FBS [Ml- In general, in these 2D lattices, the 
Mermin-Wagner theorem states that no phase transitions 
occur at finite temperatures [11. A layered structure ex¬ 
hibits specific features in the DOS as a remnant of 2D lat¬ 
tices, even though the system itself is three-dimensional 
(3D) [ll|. Here we focus on the 3D layered Lieb lattice 
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Bill, which is a test-bed for systematic investigations 
of the effects of various singularities on phase transitions. 

In this paper, we show that the magnetic transition 
temperature Tc is clearly enhanced by the FBS and VHS 
in a 3D layered Lieb lattice at half-filling using the dy¬ 
namical mean-field theory (DMFT). We determine the 
phase diagrams with several anisotropic hoppings be¬ 
tween the inter- and intralayer directions. It is shown 
that this anisotropy is a practical parameter for con¬ 
trolling TcjW for both weakly and strongly interacting 
regions. We also propose an analytical framework for 
dealing with the singularity effects with a single param¬ 
eter called a singularity weight a. We demonstrate that 
Tc oc aU for FBS and Tc/W oc for VHS. 

These forms clearly provide a unified picture of phase 
transitions dominated by singular DOSs, where a large 
singularity-weight enhances Tc- This is in stark contrast 
to phase transitions in nonsingular systems. We also 
demonstrate that, for both singularities, a conventional 
form can universally be reproduced for a —>■ 0. 

Our approach also phenomenologically explains that the 
FBS induces an anomalous behavior of thermodynamic 
quantities even in the paramagnetic state above Te¬ 


ll. MODEL AND METHOD 

Our model is described by the Hubbard Hamiltonian 
on a 3D layered Lieb lattice [see Fig. □(a)]: 

'H = —txy c\^^Cljcr—tz c|j^C/qcr-|-t/ 

l{i,j)<7 li 

where is the annihilation (creation) operator 

of an atom with spin a at site i on the /-th layer, and 
n-iitT = The subscript denotes the 

summation over the nearest neighbor sites in the xy plane 
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FIG. 1. (Color online) (a) 3D layered Lieb lattice. Solid 
lines represent the unit cell for our calculations, (b) Transi¬ 
tion temperature Tc vs. interaction U for tz = 0.1, 0.5, and 
1.0. The thick dotted lines are analytical forms: oc aU, 

oc and 2)7^® oc . See text for their full 

forms, (c) Tc vs. iz for U = 0.2. 


(z direction). We impose periodic boundary conditions 
for all X, y, z directions. 

To investigate the magnetic properties of this model 
at half-filling, we use the DMFT approach [2l|, with 
a six-site unit cell as shown in Fig.[T](a). The bipar¬ 
tite structure allows us to focus on the antiferromag¬ 
netic ordering. We employ the numerical renormaliza¬ 
tion group method (NRG) [1^ [IJ to solve the effective 
impurity problem. NRG is applicable to energy scales 
ranging from the ground state to finite temperatures Hi- 
Our approach (DMFT-I-NRG) succeeded in stud yin g 
the ground state properties of the present system 
suggesting the validity of the application to study finite 
temperature properties. The numerical procedures are 
detailed in Ref. [2l|. 

Here, we choose bandwidth W(= 4\/2txy + 4tz) as the 
unit of energy and use the notation X = X/W to de¬ 
scribe dimensionless parameters. As an exception, we 
use tz = tzjtxy to characterize the anisotropy of in¬ 
tralayer (cry-plane) and interlayer (^-direction) hoppings. 
The DOS p{ijj) is rescaled as p(w) = IFp(w). We calcu¬ 
late thermodynamic quantities, i.e., magnetization rria = 
(uaf — nai)l2 and double occupancy Da = (nafUai}, 
where a = -H,A, B. Note that the lattice geometry 
and the half-filling condition result in symmetry e.g. 
ruA = ms and Da = Db, and so on. 


III. PHASE DIAGRAM 

We first overview the characteristic magnetism of the 
present system based on the phase diagram. Figurc[T](b) 


iS 

<D 

la 




FIG. 2. (Color online) (a)-(c) Average noninteracting DOSs 
for tz = 0.1, 0.5,1.0 and (d) average magnetizations vs. 
U for tz = 0.1, 0.5, and 1.0 at zero temperature. Insets of (a)- 
(c) represent the simple DOSs obtained by extracting essential 
structures of pa.ve{to ~ 0) with normalized condition a -|- fo = 
1: p™®(cj) = aSioj) + bp™(w), P™(w) = 0{l/2 - |aj|), and 
p^*^®(cj) = aln(l/2|cj|)0(l/2 — |ca|) -|- bp'^"‘(a>), respectively. 
The inset of (d) is a closeup around t7 = 0 of for all 

three iz with analytical results (thick dotted lines): m™® ~ 
a/2 + {abU /2) ln(2/aU) with a = 0.14, /U, and 

CK with a — I/Stt. See text for the full 

form of m™® and m^'^®. 


shows Tc as a function of U for iz = 0.1,0.5, and 1.0. 
Below Tc, the antiferromagnetic insulating states appear, 
while above T)., non-magnetic (metallic or Mott insulat¬ 
ing) states appear. From nonmonotonic T). curves, we 
can see that crossovers from the band picture to the 
Heisenberg (local) picture of magnetic transitions occur 
at around U ^ 0.8-1.0 for all G. 

We next discuss how anisotropy G affects Tc- For 
clarity, we provide a change in Tc in a weakly inter¬ 
acting region U = 0.2 in Fig. [11(c). For [/< 1, Tc is 
strongly enhanced for small G(= 0.1). Surprisingly, Tc 
for tz = 0.1 shows specific behavior Tc oc U, which is 
qualitatively distinct from the well-known conventional 
weak-interacting behavior Tc oc exp(—l/G) (as found in 
those for tz = 0.5). We find that T). also slightly increases 
for G —>■ 1, which is the result of another distinct behav¬ 
ior of Tc oc exp(—l/v^) [I|. These qualitative changes in 
the magnetism will be discussed in detail below with an 
analytical approach [see thick dotted lines in Fig.|T](b)]. 

For the strongly interacting region C/ > 1 in Fig.[I](b), 
we find that Tc is enhanced as tz decreases. In contrast 
to the above, this enhancement can be understood quan¬ 
titatively: Tc can be scaled by the effective Heisenberg 
parameter Jnei- Interestingly, in this localized-spin pic¬ 
ture region, the difference in the number of adjacent sites 
plays an important role, as discussed in Sec.fVl 
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IV. WEAKLY INTERACTING REGION 


The magnetism in this region will be characterized by 
the band structure, and, in particular, the DOS at Fermi 
energy w = 0 is an important quantity. In Fig.[2](a)-(c), 
we thus provide yOave(w)[= Y.a=H,A,B Pa{^)/^] U = 0. 
As is well known, the 2D Lieb lattice = 0) has a flat 
band 0 , and therefore DOS has the delta function sin¬ 
gularity at the Fermi energy (not shown). For the present 
systems with a finite 0 ( 7 ^ 0), the flat band becomes dis¬ 
persive with a width of From Fig.[5](a), we can see 
that a large DOS still appears for a small This remain¬ 
ing flat band structure is naively regarded as the origin of 
the specific behavior Tc oc U. As shown in Fig.[2](b), the 
DOS at w = 0 decreases with increasing and then the 
FBS disappears, resulting in the conventional Tc oc 
for tz ^ 0.5. Figure[2](c) shows another interesting fea¬ 
ture of the DOS for = 1: a logarithmic singularity at 
w = 0 like a 2D VHS, leading to Tc oc [l|. 

These results suggest that, generally, a singularity of 
DOS changes the functional forms of Tc and then drasti¬ 
cally enhances Tc- Some previous studies have discussed 
logarithmic VHS effects in the 2D square lattice P, [ 2 ^ 
and also investigated the FBS effects P, [ 2 ^. However, 
to the best of our knowledge, a general analytical formal¬ 
ism for dealing with singularity effects has not yet been 
established. 

Here, we propose a general approach for revealing the 
singularity effects, which can explain the U dependence of 
Tc in Fig.[T](b). We consider the mean-field gap equation 


1 

U 



Pave (^ ) 


tanh 


f \ 

^ 2T )’ 


( 1 ) 


where pave(w) and Aave are the average DOS and spec¬ 
tral gap, respectively, with respect to sites A, B, and H. 
Gap Aave can be rewritten as Aave = Unlive, where the 
average magnetization ruave = Sa=A b 1 ™“!/^- 
simplify the multiband structure and the specific lat¬ 
tice structure, which leads to the above site-averaged 
gap equation. The average DOS can be set to a simple 
sum of the singular and the nonsingular parts: pave(A) = 
ap^{Cd) + bp^^{Cd) [see insets of Fig. [21(a)- (c)]. Here, we in¬ 
troduce the specific parameter a ( 6 ) defined as a weight of 
the singular (nonsingular) DOS with normalization con¬ 
dition a + b = 1. We simply set as the uniform DOS 
p™‘(u)) = 0 ( 1/2 — jwl), where 0 (a;) is a step function. 

Here, we should comment that our approach provides 
a general extension of the conventional forms of Tc and A 
[oc divergent p{Ep) cannot parameterize 

Tc and A any more, and instead of this, the singularity 
weight a determines these physical quantities. Note that, 
generally, any singularities of p(w) should disappear in 
an integral and a weight a is always definable: Namely, 
f p(u!)du!(= a + b) = 1 even though € M; p(w) = 00 . 
Thus, our approach is applicable to any singularities on 
any lattice geometry. 


In what follows, we show that the above simplified ap¬ 
proach with a parameterized singularity can capture the 
essence of the magnetic transition for D <C 1. We start 
with a general discussion with a of any value, which qual¬ 
itatively explains the U dependence of Tc- After that, 
with a given a, we quantitatively compare the analyti¬ 
cal form with the numerical results. We also show that 
the introduction of singularity weight a allows us to phe¬ 
nomenologically understand the anomalous behavior of 
some thermodynamic quantities. 

We first discuss the linear-D behavior of f/ shown in 
Fig.dKb). We here consider a DOS with the FBS given 
by p™®(w) = a5{Q) -I- bp™'{Cd) shown in the inset of 
Fig.|2](a). By solving Eq. (|T|) with Aave = 0, we obtain 
the transition temperature (see Appendix [^ : 

f™ = «/4tw(£^e‘/»). (2) 

where 7 is the Euler constant and ^{x) is the Lam¬ 
bert function defined as x = W(x)e^(®\ Eor a = 0, 
given yV(x) ~ x for x ^ 0, Eq. ([2]) reproduces the 
conventional behavior without the singularity: = 

j-K. Eor f/ —>■ 0, except for a = 0, the diver¬ 
gent argument requires another asymptotic prop¬ 

erty W(x) ^ Inx —In(lnx) for x —>■ -|-oo. We thus obtain 
^ aUj^T (a6C7^/4) ln(4e'’'/a7rC7) for C7 ~ 0. This 
explains 7/ oc atJ shown in Fig.[T](b) and the strong en¬ 
hancement of Tc in Fig. HKc). Importantly, this asymp¬ 
totic behavior with a divergent term indicates that, even 
if a singularity weight a is very small, the FBS changes 
the nature of the transition at around C/ ~ 0 [s^. 

We next describe the enhancement of Tc due to 
the VHS. Here we consider the DOS p)^^®(w) = 

aln(l/2|a;|)0(l/2—|u;|)-|-5p''“(a;) [see Fig.|2](c)], and then 

we obtain T/y^® ^ j-K (see Appendix IbI) . 

For a = 0, T/^® also reproduces T/f®. For C7 —7 0, 

T[y^® oc ^ which causes the higher transition 

temperatures shown in Fig.[T](b) and (c). The expo¬ 
nential decay for C7 —7 0 suggests that the enhancement 
caused by the VHS is much weaker than that of the FBS 
(see Appendix (Cl) . 

We further provide the analytical form of mave(= 
Aave/G) for any a. By solving Eq. (|T]) at T = 
0 with we obtain written 

as = al2bUW{§^e^/'^^) and = 

^ j -I I 2a _ TT^ . — 

ga aV 6i>^ /[/ (see Appendix and |B|). Both 

m™®y^Q and t=o reproduce the non-singular limit 
^ave T=o = jU for a = 0. Here, we should note 

that m™® ^ a/2-\- {abU 12) ln(2/aC/) for 17 —0, and the 
constant term a/2 explains the specific feature of the flat 
band magnetism: to™® 2 ’^q shows a jump at an infinites¬ 
imal 17(~ -1-0) [2lj. 

To show the validity of the above qualitative discus¬ 
sions, we quantitatively compare the DMFT calculations 
with the analytical results. Figurc|2](d) shows the aver¬ 
age magnetizations Wave at T = 0 calculated with the 
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DMFT. The inset in Fig.[2](d) shows that the analytical 
forms of niave T=o a = I/Stt agree well 

with the DMFT calculations. Later we will discuss the 
nonmonotonic behavior of magnetization for tz =0.1 in 
Fig. [2](d). Here, we should note that the DMFT does not 
use the simplified average DOS, suggesting the validity of 
our simplification with the extraction of the singularity. 

Here, we should discuss what determines the singu¬ 
larity weight a. For VHS systems, a can be obtained 
from the series expansion of pave(w) at around w ^ 0: 
the present system with tz = 1 has a ^ I/Stt. For FBS 
systems, the averaging assumption gives a as follows: A 
simple example is the 2D Lieb lattice {iz = 0) with a 
of 1/3, where one of the three bands is the flat band lo¬ 
cated at w = 0 0. On the other hand, for the 3D Lieb 
lattices {iz ^ 0), a is zero because all the bands are dis¬ 
persive. However, as shown in Fig.[2](a), there is a very 
narrow band at around the Fermi energy for small tz- 
Within the framework of the static mean-field approxi¬ 
mation [Eq. (P)] , the narrow band can be regarded as a 
flat band when U becomes larger than the bandwidth 
of Atz (see Appendix |D|) . To effectively explain such phe¬ 
nomena, we redefine the singularity weight a as a function 
of the other parameters: a —>■ a(f7, T, iz). 

Figure[2](d) shows that, for D > 0, ruave for iz = 0.1 
rapidly increases from zero without a jump, which can be 
phenomenologically explained by an increase in a from 
zero to a finite value as discussed above. Then, for 
[/ > 0.1, Wave increases linearly, which can be effectively 
explained by the analytical form of with a ~ 0.14 

as shown in the inset of the figure. These findings are con¬ 
sistent with Tc behavior at around U = 0 and Tc oc all 
shown in Fig.[T](b). The behavior of m-ave and Tc is char¬ 
acteristic of a flat-band magnetism in the 3D layered Lieb 
lattice [llj. 

Our DMFT calculations in Fig. P (b) and (c) clearly 
elucidate that, for small < 0.1, Tc is strongly enhanced 
by the effective FBS resulting from a narrow dispersive 
band with the interaction effects. We should stress that 
this effective FBS can be seen in various systems, such as 
the multiorbital systems with different bandwidths [s^, 
and may greatly enhance Tc of magnetism and also su¬ 
perconductivity in these systems. 

We further demonstrate that the introduction of 
a{U,T,tz) effectively explains the behavior of the ther¬ 
modynamic quantities. In Fig.P we show the average 
double occupancy Dave = J2a=H.A,B for iz = 

0.1,0.5, and 1.0 with U = 0.2. A kink in Dave clearly 
shows the transition between magnetic insulating and 
non-magnetic metallic phases. At low temperatures, 
Dave increases with increasing T; 9Dave/5T > 0 for all tz 
[33j. At higher temperatures, in the metallic region, we 
find 9Dave/9T > 0 for iz = 0.1, whereas 9Dave/9T < 0 
for tz = 0.5 and 1.0. The inset shows that, for smaller 
U = 0.1, the metallic region shows dD^ve/dT < 0 for 
all iz- In fact, dDa^g/dT < 0 can be understood from 
the usual Fermi liquid behavior, while dDa^g/dT > 0 is 
unusual as discussed below. 



T 

FIG. 3. (Color online) Average double occupancy Dave vs. 
temperature T with U = 0.2 for iz = 0.1, 0.5, and 1.0. The 
inset shows that for U = 0.1. 
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FIG. 4. (Color online) (a) Magnetization \m.y\ (7 = H,A) 
and (b) the difference of magnetization Am = \mH\ — |mA| 
vs. T with iz = 0.1 for U = 0.2, 0.8, and 1.2. 


The thermodynamic relation provides dD/dT = 
—dS/dU, where S is entropy. The Fermi liquid 
obeys S oc M*T, where AI* is the effective mass. 
The interaction-induced mass renormalization means 
dM*/dU > 0, which leads to i9Dave/(?T < 0 [2^. The 
FBS breaks down the above scenario. The entropy is 
given as oc aln2-|-6M*T, and thus cIDave/9T oc —dajdU 
at low temperatures. We conclude that the unusual be¬ 
havior dDgye/dT > 0 at D = 0.2 and iz = 0.1 results 
from da/dU < 0, meaning that the renormalization ef¬ 
fects greatly reduce the weight of fermions in the flat 
(very narrow) band at the Fermi energy. On the other 
hand, as discussed above, until U becomes comparable 
to 4:tz, we can expect da/dU > 0, which is consistent 
with what is shown in the inset of Fig.|3 Thus, our phe¬ 
nomenological approach explains the unusual Dave be¬ 
havior caused by the FBS. This unusual behavior signals 
the onset of the specific magnetic transition with Tc oc U. 

We finally discuss m-ave near transition temperatures. 
For Dc "C 1, the gap equation © provides = 

{2V3fc/U)[ifc - f)/fc\^^^ except for a = 0, and 

<ve®T^T. = {2V2iifc/Jimu) [{fc-f)/fc ]for any 
a, where (/(a:) is the Riemann zeta function (see Ap¬ 
pendix and [B]). The magnetization is proportional to 
(Tc —T)^/^ even with the singularities, which can be con¬ 
firmed by the following DMFT calculations. Figure|3](a) 
shows m..y(7 = H,A) for U = 0.2,0.8, and 1.2 with 

= 0.1 as a function of T. 
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V. STRONGLY INTERACTING REGION 

The magnetism in this region is effectively discussed 
within the local spin picture. Thus, an important quan¬ 
tity is the site-dependent rja, namely the number of ad¬ 
jacent sites (i.e., coordination number) in the xy plane 
of the site a{= A, B, H). Note that the following discus¬ 
sion will be generally applicable to bipartite lattices with 
different r]a. 

Employing a simple mean-field approach, we can ob¬ 
tain Tc = JHei /2 and J^ei = {‘^^VAVnily+^tD/U, which 
explains the quantitative change in Tc in Fig.[T](b). Here, 
Tc obeys oc Jnei = (4-\/2-|-4t^)/(4-\/2-|-44)^(l/i7), which 
takes its minimum value at around = 1 . 0 . 

We also obtain near with the mean-field ap¬ 
proach: mA(H) = meaning that 

the site with a large rja shows a large \ma\. Namely, 
|m//| > \mA\ for 17 > I at around 27 1 which is confirmed 
by the DMFT calculations shown in Fig.|3](a). In con¬ 
trast, we find |mi/| < \mA\ at low temperatures. In this 
region, the quantum fluctuations caused by the itiner¬ 
ancy of electrons yield ? 7 Q,-dependent double occupancies 
Dh > T> 4 , and a large Da suppresses the development 
of \ma\ [m|- This causes the crossing curves tuh-T and 
niA-T seen in Fig.|4](a) for U = 1.2. Furthermore, we can 
stress that the change in the sign of Am = |m^f| — jm^l 
is a clear manifestation of the crossover of magnetism 
between band and Heisenberg pictures [see Fig. |4](b)]. 


VI. SUMMARY 


We investigated magnetism in three-dimensional lay¬ 
ered Lieb lattices and determined the phase diagrams 
using the dynamical mean-field theory. We revealed that 
the (delta-functional) flat-band and (logarithmic) Van 
Hove singularities affect phase transitions, which greatly 
increases the transition temperatures for weakly interact¬ 
ing region. We also pointed out that the effective flat- 
band singularity emerges from a dispersive flat-band as a 
consequence of the interaction effects, which can appear 
in multiband systems. For strongly interacting region, we 
characterize TcjW by the number of adjacent sites. The 
larger Heisenberg interaction for site H triggers the onset 
of magnetization. Stimulated by this, we proposed a suit¬ 
able quantity, the difference of magnetization, for clearly 
detecting the crossover from the flat-band to Heisenberg 
magnetism, which can be observed in cold atom experi¬ 
ments. 

We proposed a comprehensive approach for investigat¬ 
ing the singularity effects by introducing the singularity 
weight a. We derived the universal forms of Tc and mag¬ 
netization for both singularities, which offer a remarkable 
statement: a large singularity-weight induces enhanced 
Tc- Thus, we elucidated a common feature between the 
flat-band and Van Hove singularities, which suggests that 


the singularity weight a is a unified parameter for describ¬ 
ing the singularity-induced phase transitions. 
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Appendix A: Derivation of Tf®® 


We derive Tjby solving the gap equation Eq. ([T|) . 
Substituting p^®®(a;) = a6{ui) + b6{l/2— |a;|) and Aave = 
0 in Eq. ([T]) , we obtain 


1 a 

U^Wc 


dx6{x) 


tanha; 

X 



tanh X 
ax - 

X 


Using dxtaahx/x = hilAe^y/Tr) for y > 1, we ob¬ 
tain where 7 is the Euler con¬ 

stant. Then, with the Lambert function yV{x) defined as 
X = Wix)e^'^^\ we find ^ = W Finally, 

we obtain 

r™ = «/4i.w(£:^e‘/»), (Al) 


where we assume Tc ^ l.This form reproduces Tc oc all 
for any a and T), oc for a —>■ 0 as mentioned in Sec. 

lYl 

By solving Eq. m with the same type of calculations 
as the above, we can obtain the magnetization ruave at 
zero temperature: 


ml^e:T=o = « /^bUW , 


(A2) 


where we assume Aave = maveU <C 1. It should be noted 
that shows a finite jump with infinitesimal U, 

which is a distinct feature of the flat-band magnetism. 

Employing Eq. o, C3/I1 obtS/lIl TTlg^Ye 

near the transition temperature: 


m 


FBS 

ave,T~Tc 


2y/inTc 

U 


462T 


iTc-T 


T 

A- r. 


where we use Aave = niave 
sumption 77 "C 1 , we obtain m 
except for a = 0 . 


a7r2 -h 425C(3 )Tc 

(A3) 

[/ <C 1. With an additional as- 

2 V 3 T 

u 


FBS 

ave,T~Tc 


Tc-T 


Appendix B: Derivation of 

We next derive 7)7^®. Substituting p^^®(w) = 

[aln(l/2|w|)-b 6 ] 0(1/2 — jwj) and Aave = 0 in Eq. ([TJ, 
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we obtain 


u=-\ 


^ ln(4Tca;) tanha; ^ ^ 


^i/4Te tanha; 

ax -. 


To solve the above, we give 


In a; tanha; In^ 2/ ^ ^ 

- dx = — - \- A for y >1, 


A = ji- 


■ In' 


(D 


In^ TT 




where 71 is the Stieltjes constant. Note that we use the 
following relation to derive the above constant A: 


00 

^(_l)"+i[ln2„_ln2(n + l)] = 


n—1 


7 ^ TT^ In^ 2 + ln( 4 / 7 r) InTT 

Y ~ 2 ■ 


The gap equation reduces to —aA + 61n(e'''4/7r) — [b + 
a ln(e’^/7r)] 111 ( 4 X 1 .) + ln^(4'fc) = 1/U. Then, we obtain 
the analytical form 


yVHS ^ gl+l-lVl+W + '^/TT, (Bl) 

where C = 2A + ln^(e'’' 4 / 7 r) = 7 ^ — 7 r 2/4 + 2 In^ 2 + 271 . 
The above equation holds when Tl ;< 1. 

We can obtain the magnetization Wave at zero temper¬ 
ature, by solving Eq. (ED similarly: 


m 


VHS 

ave 




(B2) 


where we assume Aave = "maveU <g; 1 . 

Using Eq. ED, we can obtain magnetization m-ave near 
the transition temperature: 


m 


VBS 

ave,T~Tc 


2^/2-kT^ 

VlU 


6 -I- 07 — aln(7rXc) 


iTr-T 


{a + 6)C(3) + aC'(3) - aC(3) ln( 7 rXc/ 26 / 7 ) V 


(B3) 


where we assume Aave = ninveU <C 1. With an ad¬ 
ditional assumption Tc ^ 1 , we obtain ^ 


2V2-kTc / Tc-T 

^ 7 C( 3 )C 7 V 


for any a. 


Appendix C: Universal ratio 

By using Eqs. (IATD-(Ib 1). we can discuss the univer¬ 
sal ratio R = 2Aave,T=o/Tc, which allows us to quantita¬ 
tively evaluate the effects of the singularities as discussed 
below. It is known that, without the singularity, R stays 
constant at 27re“'’' ^ 3.53, meaning that At=o and 
obey a similar form. Generally speaking, the singularity 
yields R depending on U and also a. 


With the FBS, the universal ratio is given by 

We find that ranges from 27re“'^ to 4 depending on 

U and a. For U <C 1, we obtain = 27re“''' (a = 0) 

and = 4 (a 0). This jump of i?^®® results from 
the divergent term in the arguments of Lambert 

functions. On the other hand, for the VHS, we obtain 


R 


iVHS 


= 27rR t/' 




-- 1 / 


l+T 


For U 1, we find i?^^® = 27re for any a. 

These findings indicate that the effects of the FBS 
(VBS) on the magnetism are strong (weak). We should 
note that, for t7 <C 1, the universal ratio depends only 
on the types of singularities, and does not depend on the 
weight a. Thus, we can stress that the universal ratio 
for (7 <C 1 characterizes the strength of the singularity 
effects. 


Appendix D: Flat band singularity caused by a 
narrow but finite bandwidth 


In this section, we explain why the narrow bands can 
be regarded as the origin of the FBS within the frame¬ 
work of the gap equation. We now consider a two- 
uniform-band (TUB) model, whose DOS is described as 
p™^(d;) = (a/IUnar)6'(H7iar/2-|a)|)-|-66>(I/2-|w|), where 
lUnar(= Hnar/lU < I) is the bandwidth of the narrower 
band Wnar scaled by that of the wider band lU. This 
DOS is naturally regarded as the simplified DOS of the 
3D layered Lieb lattice shown in Figl2](a) for U = 0.1 . 

As an example, we consider the gap equation for Tc, 
which is given by 


1 

U 


a tanhx 

W^nar Jo ^ X 



tanh X 
dx - 

X 


When both lUnar ^ 47), and 1 > 47), are satisfied, given 
dxtdxAix/x = ln(4e'’'j//7r) for y > 1, the gap equation 
produces: 


_ Wnar 

j,TUB ^ -Wna')'’""" c” C(a + bw„ar) . 

TT 

This form is qualitatively equivalent to the non-singular 
form 7)f® oc q-'^/p{Ef)u^ where p{Ep) = {a + 
WnarJ/lUnar. _ 

If Wnar ^ 4Xc and 1 > 4Xc are satisfied, by using an¬ 
other relation dxtanhx/x = y for y < 1, we obtain 

This gap equation is the same 
as that for the FBS as shown above in Appendix El 
which leads to Xj®® ~ alljA. Thus, when the interac¬ 
tion strength becomes greater than the narrower band¬ 
width, U > Wnar, we find the linear-C/ behavior of Tc. We 
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should note that here we use a relation U oc 47^. Con¬ 
sequently, we can conclude that the very narrow band 
can be regarded as the origin of the (remaining) FBS as 
mentioned in Sec. IlYl 

Here, we should comment that the above discussion 
is based on the static mean-field approximation. Thus, 
we miss some effects caused by the interactions, such as 
the mass renormalization and the creation of the Hub¬ 
bard bands. Simply put, the above gap equation can¬ 
not be used for the strongly interacting (Heisenberg) re¬ 
gion, and thus, this approach may overestimate the re¬ 


gion in which Tc (x U behavior can be obtained. Note 
that the DMFT calculations properly capture the effects 
mentioned above, and the DMFT calculations confirm 
that Tc X U behavior is obtained for small iz{^ 0.1). 

By extending the simple gap equation, we can discuss 
the above phenomena effectively. Here, we redefine a as a 
function of other parameters: a —> a(C/, T, t^). For exam¬ 
ple, we can expect that a the quasi-particle weight of the 
fermions in the narrow band will decrease owing to the 
renormalization effects. This phenomenological approach 
can explain the characteristic behavior of thermodynamic 
quantities as shown in Fig. [3] in the main text. 
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